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This article reports on a comprehensive study of the two-photon absorption (2PA) properties of six novel 
push-pull octupolar triarylamine compounds as a function of the nature of the electron-withdrawing 
groups. These compounds present an octupolar structure consisting of a triarylamine core bearing two 
3,3'-bis(trifluoromethyl)phenyl arms and a third group with varying electron-withdrawing strength (H < 
CN < CHO < N02 < Cyet < Vin). The 2PA cross-sections, measured by using the femtosecond 
open-aperture Z-scan technique, showed significant enhancement from 45 up to 125 GM for the lowest 
energy band and from 95 up to 270 GM for the highest energy band. The results were elucidated based on the 
large changes in the transition and permanent dipole moments and in terms of (i) EWG strength, (ii) degree 
of donor-acceptor charge transfer and (iii) electronic coupling between the arms. The 2PA results were 
eventually supported and confronted with theoretical DFT calculations of the two-photon transition 
oscillator strengths. 

Comprehensive studies of the two-photon absorption (2PA) cross-section on a wide spectral range are of 
fundamental importance to understand, at the molecular level, the relation between the molecular struc- 
ture and electronic and optical properties of novel organic materials. Due to its remarkable characteristics, 
such as the advantageous excitation in the near-IR region and the high spatial localization of excitation, the 2PA 
has been used in several kinds of applications including 3D optical data storage 1,2 , fluorescence imaging 3 5 , all- 
optical switching 6,7 , optical power limiting 810 , photodynamic therapy 1113 , frequency upconverted lasing 14 , coher- 
ent control of molecular systems 15,16 , microfabrication via photopolymerization 17,18 , micromachining 19,20 , to name 
a few. 

In the last three decades, different types of molecules have been proposed as alternative molecular designs to 
improve nonlinear hyperpolarizabilities and decrease the irradiance threshold necessary to generate a specific 
nonlinear optical effect. Three main designs have been largely studied: the dipolar (A-7i-D) 21 ~ 23 , quadrupolar (A- 
7t-D-7i-A) 24 ~ 26 and octupolar structures ((D-7t-A)3) 27 30 , where D and A refer to the electron-donating and 
electron-withdrawing groups respectively, linked through a bridge of 7t- conjugated bonds. On the other, some 
authors also followed the oligomeric strategy and explored the possibility of designing organic molecules 
exhibiting at the same time a strong 2PA and a good transparency 31 . Recently, great efforts have been focused 
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on the synthesis of octupolar structures since they have demon- 
strated strong electronic coupling between the individual arms, giv- 
ing rise to prodigious hyperpolarizabilities 27 . Moreover, as it has 
been shown by some of us 32,33 , even for octupolar molecules whose 
permanent dipole moment vanishes in the gas phase, aggregation 
can still take place due to the thermal motions of atoms and strong 
intermolecular interactions. 

Among the octupolar structures, push-pull triarylamine com- 
pounds stand out due to their good electrondonating and -transport- 
ing capability. Thereby, their noteworthy electronic and optical 
properties have been harnessed in different applications as active 
layers for two-photon holographic rewritable micropatterning 34 , 
fluorescence photoswitching 35 , or dye-sensitized solar cells 36 . 
However, for most of the triarylamine compounds reported so far, 
the overlap of the n molecular orbitals (MOs) carrying the nitrogen 
lone pair with the n MOs on the phenyl rings is relatively poor due to 
the large dihedral angle between the planes containing the phenyl 
ring on one side and the nitrogen-bonded carbon atoms on the other 
side 37 , hence the 2PA performances were severely impaired. We want 
here to investigate the influence of the introduction in the positions 3 
and 5 of the phenyl rings of electron-withdrawing groups such as the 
trifluoromefhyl (CF3) one, and the distinct strength of the electron- 
withdrawing groups (H < CN < CHO < N02 < Cyet < Vin) upon 
the 2PA property of triarylamine derivatives. Furthermore, the CF 3 
group imparts organic molecules with hydrophobic properties which 
exert positive influence on their use as biologically active com- 
pounds 38 . In this context, we reported here a comprehensive experi- 
mental and theoretical study on the 1PA and 2PA properties of six 
push-pull triarylamine molecules containing trifluoromethyl (CF 3 ) 
and distinct electron- withdrawing groups (EWG). Since these mole- 
cules present octupolar structures with phenyl bridges and strong 
electron-donating (triarylamine) and electron- accepting (trifluoro- 
methyl) groups, large changes in dipole moment upon excitation are 
expected. Presumably, the presence of the CF 3 may result in 
enhanced charge transfer from the amino core toward the CF 3 - 
biphenyl groups in one of the excited states. For the 2PA measure- 
ments, we employed the open-aperture Z-scan technique using an 
optical parametric amplified pumped by a femtosecond laser system 
operating at a low repetition rate (1 KHz). Two and three-energy 
level approximations helped us establish a quantitative relationship 

CF 3 




EWG 

EWG: H, CN, CHO. N0 2 



between the 2PA and the molecular properties of the chromophores 
in solution. This relationship is established by comparing the values 
of the transition dipole moments |i and dipole moment changes Au 
obtained from 2PA and solvatochromic shift measurements. The 
results were also compared with theoretical predictions of the one- 
and two-photon transitions strengths, transition dipole moments 
and dipole moment changes upon excitation, provided by response 
function calculations within the density functional theory (DFT) 
framework 39,40 . Linear optical properties such as molar absorptivity, 
one-photon induced fluorescence, transition and permanent dipole 
moments, fluorescence quantum yields and lifetimes also are 
reported. 

Results 

Synthesis of the push-pull triarylamine compounds. The synthesis 
of the targeted compounds proceeded via a Suzuki aromatic coupling 
reaction between the commercially available 3,5-bis(trifluoromethyl) 
phenylboronic acid and various bis(4-bromophenyl)phenylamino 
derivatives differing by the nature of their electron-withdrawing 
group (Figure 1). In this way, four derivatives CF 3 X where X = H, 
CN, CHO and N0 2 could be obtained as white to bright yellow solids. 
The last two derivatives CF 3 Cyet and CF 3 vin were issued from an 
additional Knoevenagel reaction between CF 3 CHO and ethyl cyano- 
acetate and malononitrile respectively, assisted by ammonium 
acetate as an acidic catalyst in a solvent mixture of pyridine:acetic 
acid (10 : 1). The compounds appear as bright orange-yellow solids. 
All six compounds were carefully purified by chromatography 
column before studies. 

One-photon absorption and emission, oxidation potentials. The 

ground-state absorption and emission spectra of the push-pull 
triarylamine derivatives CF 3 EWG are illustrated in Fig. 2a and 2b, 
respectively. All molecules present an intense absorption band in the 
UV range around 350-370 nm. A second absorption band appears at 
lower energy in the visible range for compounds CF 3 N0 2 , CF 3 Cyet 
and CF 3 vin. Interestingly, compared to their homologous derivatives 
tBuEWG, substituted by electron-donating fert-butyl groups the 
para-phenyl positionfhe CF 3 EWG series, substituted by two electron- 
withdrawing trifluoromethyl groups in the meta positions, exhibit 
bathochromic and hypsochromic shifts of their UV and visible 




EWG: H CN CHO N0 2 CHC(C0 2 Et)(CN) CHC(CN) 2 

CF 3 H CF 3 CN CF 3 CHO CF 3 N0 2 CF 3 Cyet CF 3 Vin 


increasing charge transfer in the octupolar CF 3 EWG series 




Figure 1 | Synthetic pathway toward the octupolar push-pull triarylamine compounds CF 3 EWG: i- 3,5-bis(trifluoromethyl)phenylboronic acid, 
Pd(PPh 3 ) 4 , Na 2 C0 3 2 mol.L -1 , toluene, 80°C, 8 h; ii- ethyl cyanoacetate or malononitrile, pyridine:acetic acid (10 : 1), NH 4 OAc, RT, 6 h. 
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Figure 2 | (a) Ground-state absorption and (b) fluorescence spectra of push-pull triarylamine derivatives containing CF 3 in toluene solutions. 



absorption bands respectively. However, their molar absorption 
coefficients e are similar to those obtained for the tBuEWG series 41,42 . 

As for the absorption spectra, only one fluorescence band is 
observed in the visible region. Interestingly, two kinds of photophy- 
sical behaviors can be distinguished. The first one relates to com- 
pounds CF 3 N0 2 , CF 3 Cyet and CF 3 vin. All three display an 
emission band centered around 520 nm and a large Stokes shift ran- 
ging from 3900 to 6300 cnT 1 . As already discussed for the tBuEWG 
series 43 , femtosecond transient absorption spectroscopy measure- 
ments showed that the compounds underwent a significant nuclear 
rearrangement in the excited state after a charge transfer occurred 
from the triphenylamino core onto the EWG unit. The second beha- 
vior is described by compounds CF 3 H, CF 3 CN and CF 3 CHO which 
exhibit an emission band centered around 430 nm and similar Stokes 
shifts ranging from 4300 to 5100 cnT 1 . Compared to the other two 
compounds, the emission band of CF 3 CHO lies at lower energy in 
accordance with the more red-shifted wavelength of the correspond- 
ing maximum absorption. Finally, the emission of CF 3 H in toluene 
located at 428 nm instead of 383 nm for the tBuH homologue shows 
a clear participation of the trifluoromethyl groups in the radiative 
excited state and presumably in the Franck-Condon state also. 

The electronic influence of the CF 3 groups is clearly apparent in 
the oxidation of the triphenylamino core into a triphenylammonium 
cation as featured by the first redox potential (Table 1). All triaryla- 
mine derivatives CF 3 EWG examined in this study displayed a 



quasi- reversible process characterized by an averaged redox poten- 
tial E° ranging from 0.63 to 0.87 V (vs Fc + /Fc) depending on the 
EWG substituents). Compared to the series of electron-donating 
compounds tBuEWG, the E° values rise by 0.17 to 0.2 V. This con- 
siderable increase, being much more significant for CF 3 H, is to be 
ascribed to the electron-deficient character of the CF 3 groups, thus 
decreasing considerably the electron density located on the triphe- 
nylamino core in the ground state. No simple correlation between the 
electrochemical measurements and the optical gaps could however 
be traced, which is due to the high sensitivity of the excited state on 
the solvent polarity as it will be discussed later. 

Character of electronic transitions and solvatochromic measure- 
ments. In order to investigate the nature of the transitions observed 
along the 1PA spectra, we performed the analysis of respon- 
se vectors to determine the largest one-electron orbital excitations. 
As it can be noted in Fig. 3, the lowest-energy transition for the H and 
CN derivatives is a n — >tc* transition involving the two trifluoro- 
methyl-phenyl arms (acceptor groups). For all other compounds, 
the lowest-energy transition is an ICT transition involving the triary- 
lamine core and the respective EWG. The reason of this ordering 
difference stems in the low strength of the CN substituent, which 
therefore does not promote the inversion between the it* and ICT 
excited states as observed for all the other substituents. The influence 
of the EWG strength on this inversion is also evidenced by the 
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Figure 3 | KS molecular orbitals involved in the lowest-energy electronic transitions of the investigated molecules. 



molecular orbital involved in the excitations of the CHO substi- 
tuent. Because of its moderate EWG strength, the charge transfer 
from the triarylamine core to the EWG in the case of CHO substi- 
tuent is not as pronounced as for the N0 2 , Cyet and Vin substituents. 

To shed more light on the character of electronic transitions, we 
performed solvatochromic Stokes shift measurements to find the 
difference between the permanent dipole moments in the excited 



and ground states, A/i 01 . 
Mataga equation 44 : 



. ^cci2 3nhc dv 



Vol 



(1) 



where u = v A —u em is the difference between the wavenumbers of the 
maximum absorption and fluorescence emission (in cm— 1), 
F(n,f) = 2(<i;-l)/(2£ + l)-2(n 2 -l)/(2ra 2 + l) is the Onsager polarity 
function with Ej being the dielectric constant of the solvent, Vol is the 
volume of the effective spherical cavity occupied by the molecule 
inside the dielectric medium. Various solvents or mixtures of sol- 
vents (cyclohexane, hexane, toluene, hexane/chloroform (50/50%), 
hexane/dichloromethane (50/50%), chloroform, dichloromethane, 
acetone and ethanol) with increasing polarity were investigated for 
all compounds. 

A pronounced increase of the Stokes shift was found with increas- 
ing solvent polarity and reveals clearly the polar character of the 



radiative excited state (Si) (Fig. 4). These results indicate that the 
absolute values of the excited state dipole moments p, u are much 
higher than those of the ground state dipole moments Jl Q0 . 



ovh _, dv. 
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Using these data, the hydrodynamics volume estimated from the 

IEF-PCM/B3LYP/6311G(d,p) calculations (see Table 1) and the 

Eq. 3, we found |A/J^| = 10.6D, |A$jf|=9.2D, |Aju£ 

and 



9.5D, 



I A -NO, 



= 12. 5D, 



A/i 0 C f =14.7D 



|A/1oTI = 



Ll 01 I - >■'■■■»->■< <-»/-'01 — lt.iis aim |*-V*01 

14. 8D. These results show that all molecules present pronounced 
A(i 01 values (between 9 to 15 Debye) indicating a high degree of 
effective charge transfer (>34%), which can be obtained from the 
relation r e ff = \ Ay\/e, where e is the electron charge 45 . Table 1 reports 
the degree of charge transfer observed for all molecules. 

As one can note, the degree of ICT (jeff / r)> where r is the distance 
between the core and the first EWG atom obtained from the DFT 
calculations, for the molecules CF 3 CN and CF 3 CHO is considerably 
smaller than that of CF 3 H, even though the cyano and aldehyde 
substituents have a higher EWG strength. Another important out- 
come that we can infer from these results is the increase in the degree 
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Figure 4 | (a) Solvatochromic Stokes shift (u) measurements obtained as a function of the Onsager polarity function (F(n, {)) for all molecules in a series 
of different solvents (cyclohexane, hexane, toluene, hexane/chloroform, chloroform, dichloromethane, acetone and ethanol). (b) Solvatochromic 
measurements for CF 3 Cyet molecule illustrating the increase of the Stokes shift with the increase of the solvent polarity for this class of molecules. 



of ICT caused by higher splitting of the lowest energy 1PA band. The 
splitting in quasi-degenerate excited states is a result of the interac- 
tions between the arms of the chromophores 27 . Whereas no such 
splitting could be experimentally distinguished for CF 3 CN and 
CF 3 CHO that have the lowest values of ICT, clear interactions could 
be noticed for the other four compounds: a shoulder on the lowest 
energy band for CF 3 H with an energetic separation of —300 meV, 
and a considerable splitting for CF 3 N0 2 (~460 meV), CF 3 Cyet 
(-610 meV) and CF 3 Vin (~740 meV). We used the Gaussian 
decomposition method to estimate the splitting energy for CF 3 CN 
and CF 3 CHO and obtained 240 meV and 250 meV, respectively. 
From these values and based on the exciton Frenkel model 27 , we 
estimated the electronic coupling constant for all octupolar struc- 
tures as being —100 meV for CF 3 H, —80 meV for CF 3 CN and 
CF 3 CHO, -155 meV for CF 3 N0 2 , ~205 meV for CF 3 Cyet and 
—245 meV for the CF 3 Vin. As a conclusive point, our results show 
that the electronic coupling is directly proportional to the degree of 
ICTofEWG. 

Table 1 also reports the photophysical parameters of push-pull 
triarylamine derivatives containing CF 3 , such as: maximum absorp- 
tion and fluorescence wavelengths, maximum molar absorptivity, 
experimental ground-state ( | /^oo I ) an d transition dipole moment 
( I jt2 01 1 ), fluorescence quantum yield and lifetime determined in this 
work by using a series of spectroscopic techniques. All molecules 
show a significant fluorescence quantum yield with average lifetimes 
of the order of 2 ns. 

Experimental 2PA properties. Figure 5 shows the experimental 2PA 
spectra (dots) determined by performing open-aperture Z-scan 
measurements with femtosecond pulses. The experimental data 
reveal well-defined 2PA spectra with moderate cross-section values 
(from several tens to few hundreds of GM units) in the visible and 
infrared regions depending on the electronic influence of the 
withdrawing substituents. 

Along the 2PA spectra we can distinguish, except the molecules 
with the lowest degree of ICT (CF 3 CN and CF 3 CHO), three 2PA 
allowed bands and the resonance enhancement effect 46,47 . Among the 
three 2PA bands, those two ones of lower energy are related to the 
ICT (So^Si) transition and the 71-71* (S 0 ^S 2 ) transition located on 
the trifluoromethyl-phenyl arms in agreement with the linear 
absorption spectra. This assumption is supported by the good spec- 
tral correspondence observed between the bands in the 2PA and 1PA 
spectra (see Fig. 5) and also by the results of the response function 



calculations (see theoretical results section). Slight blue shifts 
between the spectral position of 1PA and 2PA peaks may be ascribed 
to the distinct vibronic coupling between the ground- and excited- 
states in the 1PA and 2PA transitions. Although the computational 
strategies to include vibronic coupling are well established 48 , we have 
not attempted to apply them in the present study. The primary 
reason is that the experimental 2PA spectra do not exhibit a fine 
vibrational structure for the absorption bands. However, in the case 
of multi-branched structures, vibronic contributions to the two- 
photon absorption may be substantial 48 . It is important to mention 
that two-photon transitions are different from the one-photon ones 
and, moreover, very strict. However, in compounds that present a 
low molecular symmetry, so that it is not possible to define the parity 
of excited states precisely, the same electronic state may be accessed 
via both one and two-photon absorptions 49 . 

It is worth mentioning that, for all molecules, the highest energy 
2PA band has the highest cross-section values that are between 
95 GM (for CF 3 H) to 270 GM (for CF 3 Vin). This band is associated 
with electronic states that present high 2PA probability located in the 
UV region (from 260 nm to 300 nm). 

The cross-section values, at the peak of the main 2PA bands, were 
determined to be c.a.65 GM at 725 nm, 83 GM at 620 nm and 
97 GM at 520 nm for CF 3 H; 58 GM at 715 nm and 145 GM at 
560 nm for CF 3 CN; 48 GM at 720 nm and 262 GM at 580 nm for 
CF 3 CHO; 45 GM at 800 nm, 58 GM at 680 nm and 230 GM at 
585 nm for CF 3 N0 2 ; 115 GM at 850 nm, 68 GM at 700 nm and 
250 GM at 610 nm for CF 3 Cyet; 125 GM at 880 nm; 70 GM at 
700 nm and 265 GM at 605 nm for CF 3 Vin. The 2PA cross-section 
values reported here for the triarylamine molecules containing CF 3 
are at least 30% higher than those reported for their ferf-butyl ana- 
logues 41,42 . However, we cannot attribute this 2PA enhancement to 
the CF 3 groups because this difference may be associated only with 
the higher polarity of solvent used in this present report (toluene) 
when compared to the results of Refs 40,41 . (hexane). To overcome this 
doubt, we performed 2PA measurements using the Z-scan technique 
for the ferf-butyl analogues in toluene solution and, in fact, observed 
2PA cross-section values similar ones reported here to the molecules 
containing CF 3 groups. 

Analogously to what was found in the 1PA spectra of CF 3 CN and 
CF 3 CHO, no splitting between the two lowest energy 2PA transitions 
is observed along the nonlinear spectra. Additionally, the peak of 
lowest energy 2PA band is slightly red-shifted (50 meV) as com- 
pared to the corresponding 1PA peak, indicating that the S 0 — »Si 
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Figure 5 | The solid lines and dots represent, respectively, the linear and two-photon absorption spectra. 2PA cross-sections for the CF 3 H, CF 3 CN, 
CF 3 CHO and CF3NO2 molecules were multiplied by 2 for better comparison with the others molecules. 



transition dominate the nonlinear optical effect in this region. This 
result can be attributed to the low electronic coupling constant 
observed for these molecules (—80 meV). 

For the other molecules, as it was observed in the 1PA spectra, the 
2PA spectra exhibit a clear and increasing (in energy) splitting of the 
lowest energy band as a function of increase of the degree of ICT. 
Our results point out that for CF 3 H (af PA / 'of PA = 1.3) and CF 3 N0 2 
(af PA /a s 2PA = l.2), the n-n* (S 0 ^S 2 ) transition has higher 2PA 
cross-section values than that of the corresponding S 0 ^Si transition, 
whereas for CF 3 Cyet (of PA /of PA = 0.59) and CF 3 Vin (of PA /of PA = 
0.56) the opposite is found. This result can be at least partially 
explained due to the lower transition dipole moment (|/%|)> 
observed for CF 3 H and CF 3 N0 2 as compared to CF 3 Cyet and 
CF 3 Vin. In the context of the sum-over essential states approach, 
the 2PA cross-section of the lowest energy transition, which is 
proportional to the product of |/I 01 | 2 and |A/I 01 | 2 , became less 
intense. 



Since the 2PA cross-section of the S 0 ^S 2 transition for CF 3 H is 
the highest one despite a larger ICT encountered for the other octu- 
polar compounds CF 3 N0 2 , CF 3 Cyet and CF 3 Vin, one can suspect 
for CF 3 H a large electronic coupling between its bis(trifluoro- 
methyl)biphenyl arms upon 2PA excitation. The lower 2PA cross- 
sections for CF 3 N0 2 , CF 3 Cyet and CF 3 Vin molecules may be due to 
increased angular torsions of the N0 2 , Cyet and Vin substituents 
around the phenyl-N bonds, thus leading to a reduced electronic 
coupling between the arms and, as a consequence, a decrease in 
the 2PA cross-section of S 0 ^S 2 electronic transitions. 

To establish a more quantitative relationship between the effect of 
distinct EWG and the observed 2PA activity, we used the sum-over- 
essential states approach to determine the permanent dipole moment 
change between the ground and excited states and the transition 
dipole moments between excited states. Details about this approach 
can be found in "Supporting information". Figure 6 shows a 
comparison between the |A/Iq^ a | values obtained from 2PA and 
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Figure 6 | Comparison of the permanent dipole moment change (A/J 01 ) 
obtained from 2PA and solvatochromic Stokes shift measurements. 



solvatochromic shift measurements. As it can be observed in Fig. 6, 
there is a good agreement between the results determined from 2PA 
and solvatochromic shift measurements. These results show that the 
lowest-energy 2PA band for all molecules can be well-described by a 
two-level system. 

To model the S 0 ^»S 2 2PA allowed transition, we considered a 
three-level system composed of the ground state (S 0 ), one intermedi- 
ate 1PA allowed excited state (| 1)) and the 2PA allowed final excited 
state (S 2 ). Table 1 shows the results obtained using the methodology 
described in "Supporting information". The main aspect to be high- 
lighted, based on these results, is that upon going from the neutral 
substituent (H) to the strongest EWG (Vin), the transition dipole 
moment that connects the two lowest energy 2PA allowed states, 
|/2 12 |, decreases by about 2.5-5 times, due to the influence of the 
EWG on the electronic coupling between the arms of the octupolar 
structure. For CF 3 CN and CF 3 CHO, a decrease in the transition 
dipole moment |/t 12 | was expected since these molecules present 
lower electronic coupling constant as compared to CF 3 H. For the 
other molecules, the N0 2 , Cyet and Vin substituents should increase 
the angular torsion of the whole plane phenyl-EWG along the N- 
phenyl bond, which impairs electronic coupling between the arms, 
thus decreasing \p. n \- In this context, the structure-property rela- 
tionships outlined in this study are of foremost importance for the 
development of novel organic materials with remarkable nonlinear 
optical properties. 

Theoretical results. In what follows, we present the results of 
quantum- chemical calculations of the electronic structures of triary- 
lamine molecules using response functions. The data presented in 
Table 2 correspond to the results of calculations performed using 
reparametrized CAM-B3LYP functional and some additional infor- 
mation extracted from the experimental data. In general, a satisfac- 
tory agreement between the transition energies determined by the 
response function calculations and through the theoretical fitting of 
the experimental 1PA spectrum is observed. The difference found for 
the lowest energy transition S 0 ^Si is smaller than 0.15 eV for all 
molecules. In the case of the second transition reported in Table 2, 
So^S n=2i 3, the difference observed is smaller than 0.18 eV basically 
for all molecules, the only exception (0.26 eV) being the case of CN 
substituent. 

The states accessed via two-photon transitions are also allowed via 
one-photon transition as shown by the reasonable oscillator strength 
obtained for all the transitions (see Table 2). Such results confirm 
that the electric-dipole selection rules are relaxed for all the 



compounds investigated, which was expected due to the absence of 
centrosymmetric. 

By comparing the results of the gas-phase (not shown) and PCM 
calculations, the solvent effects promote a red-shift and the increase 
in the oscillator strength for all transitions, the only exception being 
the second transition of Cyet substituent that shows a small blue- 
shift. In addition, the comparison shows that the solvent effects 
promote a 2PA probability increase for all the transitions of the 
investigated substituents. The solvent effects increase the 2PA prob- 
ability by about 80-130%. To better compare the experimental data 
and theoretical results for 2PA process, we simulated the 2PA cross- 
section spectrum of the triarylamine molecules by using the excited 
state linewidths estimated by fitting the experimental linear spectra 
and the 2PA transition probabilities determined by the quadratic 
response calculations (see Table 2). The simulated and experimental 
2PA cross-section spectra are presented in Figure 7. As it can be seen, 
there is a good agreement between the simulated and experimental 
spectra both for the spectral behavior and cross-section magnitude. 

For each investigated molecule, the simulated 2PA cross-section 
spectrum reproduces the bands observed over the experimental spec- 
trum. The magnitude of the 2PA cross-section at the peak of each 
absorption band observed in the experimental spectrum and the 
related value in the simulated spectrum are gathered in Table 2. 

Discussion 

The experimental 2PA cross-sections observed along the nonlinear 
optical spectra for the CF 3 EWG molecules as a function of the nature 
of the electron-withdrawing group can be explained from the com- 
petition among the EWG strength, the degree of donor-acceptor 
charge transfer and the electronic coupling among the arms present 
in each molecule. Based on these parameters, our experimental 
results showed that the Cyet and Vin substituents favors the first 
2PA band (which has the ICT character) in detriment to their second 
band due to the poor electronic coupling between the arms and the 
strong donor-acceptor (N-^EWG) charge transfer. However, 
CF 3 Cyet and CF 3 Vin have an electronic coupling constant higher 
than those observed in the case of CF 3 H and CF 3 N0 2 which present 
2PA cross-sections higher to the second 2PA band as compared with 
the first one. Nevertheless, in solution, the octupolar CF 3 Cyet and 
CF 3 Vin molecules present a considerable increase in angular tor- 
sions of the whole plane phenyl-EWG along the N-phenyl bond that 
lead to a reduced electronic coupling between the arms, thus decreas- 
ing the 2PA cross-section of the S 0 ^S 2 electronic transition. These 
angular torsions were recently observed through femtosecond tran- 
sient absorption spectroscopy measurements for similar molecules. 
On the other hand, for CF 3 CN and CF 3 COH, the fact of the So^Si 
transition dominates the lowest energy 2PA band (at least 70%) is 
due exclusively to the low electronic coupling constant observed for 
these molecules (—80 meV). 

All CF 3 EWG molecules present the highest energy 2PA band with 
the highest cross-section value which ranges from 95 GM for CF 3 H 
to 270 GM for CF 3 Vin. This band is associated with an cluster of 
electronic transitions that present high 2PA probability located in the 
UV region (from 260 nm to 300 nm) as well as the resonance 
enhancement effect due to the little detuning between the excitation 
energy and the energy of the lowest energy one-photon allowed state. 

By comparing the magnitude of the 2PA cross-section of the two 
lowest energy transitions, the theoretical results and experimental 
data show a nice match, i.e. in the case of the H and N0 2 substituents 
the higher energy transition presents a larger cross-section. For the 
remaining substituents, the theoretical results and experimental 
data do not agree about which transition presents the higher 2PA 
cross-section. At least two factors may dramatically affect the agree- 
ment between the theoretical results and experimental data. In the 
case of the CN and CHO substituents the two-lowest energy transi- 
tions are close in energy, as evidenced by the theoretical results and 
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Table 2 Experimental data and theoretical results of the PCM/Reparametrized-CAMB3LYP calculations using the linear and quadratic 

f i- ~r •■■ ■ t r~\ • "Hi i ill li i •■■ 1 1 ■ 1 ■ ■ / rM • 1*1 i 1 

response functions. 1 ransition energies (t) in nm, oscillator strength, two-photon transition probability (r) in a.u., dipole moment changes 
(Au) in Debye and two-photon absorption cross-section (a) in GM. The molecular orbitals involved in the excitations describing the electronic 
transitions reported are also showed. AjU^H is one obtained from the 2PA cross-section data 
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the structureless band observed in their 1PA spectrum. In such a 
situation, the accuracy of the procedure to fit the experimental linear 
spectra and estimate the spectral linewidth of the related excited 
states, which is then used to compute the 2PA cross-sections, cannot 
be completely reliable. Since for the CN and CHO substituents, the 
2PA transition probability of the two lowest energy transitions are 
roughly similar, an accurate estimative of the linewidth of the excited 
states is crucial to determine which transition presents the higher 
2PA cross-section. For the Cyet and Vin substituents, on the other 
hand, the two-lowest energy transitions are separated in energy and, 
therefore, the estimated values of the spectral linewidths are more 
reliable. Also, the 2PA probabilities of the two lowest energy transi- 
tions are quite different for these substituents. In order to find a 
plausible explanation for these differences, the electronic character 
of the lowest-energy transitions should be further investigated. 

A good agreement between theoretical results and experimental 
data concerning the efficiency of the distinct EWG in increasing the 
2PA cross-section of the Sn^Sj transition is observed. For the CHO, 
N0 2 , Cyet and Vin substituents, the cross-section ordering basically 
agrees with the EWG strength ordering (CHO < N0 2 <Cyet< 
Vin). As mentioned above, in the case of CN substituent, such a 
transition does not involve its EWG. The quantitative agreement 
between theoretical and experimental 2PA cross-section is, however, 
satisfactory. 

Because the second-lowest energy transition of the CHO, N0 2 , 
Cyet and Vin substituents is a n— >7i* transition centered on the two 
trifluoromethyl-phenyl arms of these systems, the 2PA cross-section 
of such transitions is expected to be enhanced by the coupling 
between the arms. In addition, since trifluoromethyl is an EWG, in 
principle its presence also contributes to the intensity of the 2PA 
transition centered on the arms. The higher theoretical values for 
the 2PA cross-section of the ti — >7t* transition in comparison with the 
experimental values for such substituents may feature an impaired 
coupling due to the mobility of the arms in solution. Moreover, 
by comparing the present results for the CHO, N0 2 and Vin 
substituents with the results of a previous study in which similar 
systems without trifluoromethyl groups were investigated 42 , one 
can conclude that effects of the presence of the trifluoromethyl group 
are not significant for the 2PA transition. In fact, by checking the 



electron density distribution of the molecular orbitals (Figure 7) 
involved in the excitation describing the tt — >7i* transition, one can 
observe that the trifluoromethyl groups do not take part in such a 
transition. 

Methods 

Compounds. All chemical reagents and solvents were purchased from commercial 
sources (Aldrich, Acros, SDS) and used as received. Spectroscopic grade solvents 
purchased from Aldrich were used for spectroscopic measurements. Bis(4- 
bromophenyl)phenylamine 42 , 4-di(4-bromophenyl)aminobenzonitrile 42 , 4-di(4- 
bromophenyl)aminobenzaldehyde 42 and bis(4-bromophenyl)-4-nitrophenylamine 34 
were synthesized according to previous reported procedures. Details about the 
synthesis for the CF 3 H, CF 3 CN, CF 3 COH, CF 3 N0 2 , CF 3 Cyet, CF 3 Vin molecules can 
be found in "Supporting information". 

Electrochemical measurements. Electrochemical measurements were performed 
with a three-electrode cell consisted of a glassy carbon working electrode (1.5 mm 
diameter), a counter electrode made of a Pt wire coil and an Ag/Ag-t- reference 
electrode. The Ag/Ag + reference electrode was obtained by introducing a silver wire 
in the reference compartment containing silver ions (AgN0 3 , 0.01M) and 
tetrabutylammonium perchlorate (TBAP, 0.1 M) in acetonitrile. The measured redox 
potential were expressed vs ferrocene used herein as internal reference. The 
triarylamine solutions were prepared in dichloromethane solvent added with 0.1 M 
of tetrabutylammonium hexafluorophosphate (Bu 4 NPF 6 ) as a supporting electrolyte 
(Sigma-Aldrich). The solutions were deoxygenated by argon bubbling prior to each 
experiment. Cyclic voltammograms (CVs) were performed with VSP potentiostat 
controlled by EC-Lab 10.02 software. Both potentiostat and software were purchased 
from Bio-Logic Instruments-France. 

1PA and 2PA measurements. Triarylamine derivatives/toluene solutions were 
prepared with concentrations of 1 X 10 16 molecules/cm 3 and on order of 10 18 
molecules/c-m 3 for linear and nonlinear optical measurements, respectively. The 
samples were placed in 2 mm thick quartz cuvettes for optical measurements. The 
steady-state absorption and photoluminescence spectra were recorded using a 
Shimadzu UV-1800 spectrophotometer and a Perkin Elmer LS55 fluorimeter, 
respectively. 

The nonlinear optical properties were measured by using, the open aperture Z- 
scan technique, employing 120-fs laser pulses from an optical parametric amplifier 
pumped by 150-fs pulses (775 nm) from a Ti:sapphire chirped pulse amplified 
system. The Z-scan measurements were carried out with intensities ranging from 
60 to 200 GW/cin 2 (30 to 120 nj/pulse) and with beam waist size at the focus 
varying from 15 to 18 um. To ensure a Gaussian profile for the laser beam used in 
the experiments, spatial filtering is performed before the Z-scan setup. A silicon 
detector was employed to monitor the laser beam intensity in the far- field. To 
improve the signal to noise ratio, the oscillatory Z-scan method was used such 
that the samples are continuously scanned, repeating the experiment several times. 
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Figure 7 | The dots represent the experimental 2PA cross-section spectra. The solid lines correspond to the simulated 2PA cross-section spectra 
obtained using the 2PA transition probabilities determined by the response function calculations and excited states linewidths estimated by fitting the 
experimental nonlinear spectra using the sum-over essential states approach. 



Moreover, a lock-in amplifier served to integrate 1000 shots for each point of the 
Z-scan signature. 

In the open aperture Z-scan technique, 2PA cross-section was determined by 
translating the sample through the focal plane of a focused Gaussian beam, while 
transmittance changes in the far field intensity are monitored. For a 2PA process, 
the light field created an intensity dependent absorption, a = a 0 + (31; in which I 
is the laser beam intensity, oc 0 is the linear absorption coefficient, and (J is the 2PA 
coefficient. The nonlinear coefficient [i was obtained by fitting the Z-scan data. 
The 2PA cross -section, <3 2 pa> was determined from <j 2 pA~hco^/N, where fyois the 
excitation photon energy, and N is the number of chromophores per cm 3 . The 
2PA cross-section was expressed in Goppert-Mayer (GM) units, where 
lGM=lX10" 50 cm*.s.photon~\ 

Fluorescence quantum yield and time resolved measurements. Steady-state 
emission spectra were recorded using a Perkin Elmer LS55 or a Fluoromax 3 (Horiba- 
Jobin-Yvon) spectrofluorimeters respectively. Fluorescence intensity decays were 
measured by the time- correlated single-photon counting method (TCSPC) with a 
femtosecond laser excitation at 340 or 400 nm provided by a Spectra-Physics setup 
composed of a Titanium- Sapphire Tsunami laser pumped by a doubled YAG laser 
(Millennia), pumped itself by a two-laser diode array, and doubling and tripling LBO 
crystals. Light pulses were selected by optoacoustic crystals at a repetition rate of 
4 MHz. Fluorescence photons were detected at the emission maximum through a 
monochromator by means of a Hamamatsu MCP R3809U photomultiplier, 
connected to a constant-fraction discriminator. The time -to -amplitude converter was 
from Tennelec. Pulse deconvolution was performed from the time profile of the 
exciting pulse recorded under the same conditions by using a Ludox solution. The 
fluorescence data were analyzed by a nonlinear least-squares global method using the 



Globals software package developed at the Laboratory for Fluorescence Dynamics at 
the University of Illinois at Urbana- Champaign (Globals software, Irvine, CA, USA). 

Computational details. In this paper, quantum chemical calculations were employed 
to investigate the molecular and electronic structures of triarylamine molecules 
containing trifluoromethyl and determine their one-photon and two-photon 
transitions and related electronic structure parameters. For this purpose, quantum 
chemical calculations using the Kohn-Sham formulation of density functional theory 
(KS-DFT) were performed. The equilibrium molecular geometry of the molecules 
was determined in both gas-phase and solvent environment. The solvation effects on 
the equilibrium molecular geometry, the electronic structure and the one-photon and 
two-photon properties of each investigated molecule were taken into account by 
employing a Polarizable Continuum Model (PCM) and assuming the toluene solvent. 

The gas-phase geometry optimizations were performed employing the B3LYP/6- 
31lG(d,p) approach 50-52 . In the case of the geometry optimization taking into account 
the solvent environment, the PCM using the integral equation formalism variant 
(IEF-PCM) 10 was employed in combination with the B3LYP/6-31lG(d,p) approach. 
All the geometry optimizations were performed using the Gaussian 2009 package 53 . 
The IEF-PCM method was combined with the United Atom for Hartree-Fock 
(UAHF) model to build the solute cavity. In this model the van der Waals surface is 
constructed from spheres located on heavy (that is, non-hydrogen) atoms only 
(United Atom approach) . The van der Waals radius of each atom is a function of atom 
type, connectivity, total charge of the molecule and the number of attached hydrogen 
atoms. The hydrodynamic volume of each molecule in toluene was theoretically 
determined from the PCM-B3LYP/6-31lG(d,p) calculation. These theoretical values 
were used to estimate the permanent dipole moment changes under excitations based 
on experimental data. 
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Subsequently, to determine the lowest energy one- and two-photon allowed 
transitions and related molecular parameters of the investigated molecules, the res- 
ponse function calculations 39,40 within the DFT framework were performed using the 
DALTON program 54 . In this approach, the excitation energies and transition dipole 
moments {two-photon probabilities) are analytically computed as, respectively, poles 
and single residues of the linear (quadratic) response function of the molecular 
electronic density. Moreover, by computing the double residue of the quadratic 
response function the excited states dipole moments and the transition dipole 
moment between excited states were also determined. All the response function 
calculations were performed employing the CAM-B3LYP/6-3 1 + G(d) approach 55 for 
the investigated molecules in both gas-phase and solvent environment. In the case of 
the response function calculations, the PCM as implemented in DALTON program 
was employed 56 . To be consistent, for all molecules the van der Waals surface used in 
the response function calculations was defined by adopting the set of van der Waals 
radius and atomic centers determined at the geometry optimization stage using the 
Gaussian 2009 package. 

A preliminary assessment of the exchange- correlation functional revealed that the 
standard hybrid B3LYP and PBEO functional and also the standard long-range 
corrected CAM-B3LYP functional, fail to well reproduce the experimental UV-vis 
spectra of the triarylamine molecules containing trifluoromethyl investigated (for a 
recent exhaustive summary of TD-DFT benchmarks we refer to excellent 
reviews 57-59 ). 

Therefore, in this study the successful strategy of K. Okuno et al. 32 of employing 
a tuned CAM-B3LYP functional was employed. In their systematic study, K. 
Okuno et al. 32 showed that the use of the time-dependent density functional 
theory (TD-DFT) scheme with the tuned CAM-B3LYP functional is found to 
semi -quantitatively reproduce several excitation energies identified from the 
experimental UV-vis spectra of 15 closed forms of diarylethene derivatives. In 
that study, the default value of the CAM-B3LYP functional parameters, that is jj. 
- 0.33, a - 0.19, and p = 0.46, were tuned to u = 0.150, a = 0.0799 and p = 
0.9201 and the TD-DFT calculations were performed assuming the gas-phase 
condition. In the present study the values of the functional parameters were 
defined based on the comparison between the energy of the one-photon lowest- 
energy transition determined by the PCM/Response function calculation and the 
energy of the lowest-energy absorption band of the experimental UV-vis spectra 
of the smallest system investigated, that is, the H substituent. Note that in the 
present study we neglect the zero-point vibrational energy and we compare ver- 
tical excitation energy with the experimental absorption band maximum 60 . Here, 
the following values were adopted as the CAM-B3LYP functional parameters: u — 
0.150, a — 0.03 and P = 0.97. The suitability of such values to the response 
function calculations addressing the other substituents was confirmed through the 
satisfactory reproduction of the excitation energy of the lowest-energy transitions 
(one or two) identified from the experimental UV-vis spectra. 

One may point two plausible sources of discrepancies between the values adopted 
for the CAM-B3LYP functional parameters in the present study and the one by K. 
Okuno et al. The distinct class of compounds investigated is perhaps the main reason. 
However, another point that can contribute to such differences is that in the present 
study the values adopted for the CAM-B3LYP functional parameters were deter- 
mined by performing linear response function calculations in combination with a 
PCM, while in the work conducted by K. Okuno et al. the values were determined 
based on the results of gas-phase calculations. 

Finally, to confirm the donor-acceptor charge transfers involved on the main one- 
and two-photon transitions of the substitutes, the percentage contribution of the 
excitations to the transitions, and the KS molecular orbitals involved, were deter- 
mined. The KS molecular orbitals were determined by DFT calculations using the 
tuned CAM-B3LYP/6-31 + G(d) approach in combination with IEF-PCM and using 
the Gaussian 2009 package 53 . 
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